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ABSTRACT 

Kilometric-scale optical imagers seem feasible to realize by intensity interferometry, using telescopes primarily 
erected for measuring Cherenkov light induced by gamma rays. Planned arrays envision 50-100 telescopes, dis- 
tributed over some 1-4 km 2 . Although array layouts and telescope sizes will primarily be chosen for gamma-ray 
observations, also their interferometric performance may be optimized. Observations of stellar objects were nu- 
merically simulated for different array geometries, yielding signal-to- noise ratios for different Fourier components 
of the source images in the interferometric (u, i>)-plane. Simulations were made for layouts actually proposed for 
future Cherenkov telescope arrays, and for subsets with only a fraction of the telescopes. All large arrays provide 
dense sampling of the (iz, v)-plane due to the sheer number of telescopes, irrespective of their geographic orien- 
tation or stellar coordinates. However, for improved coverage of the (u, v)-plane and a wider variety of baselines 
(enabling better image reconstruction) , an exact east-west grid should be avoided for the numerous smaller tele- 
scopes, and repetitive geometric patterns avoided for the few large ones. Sparse arrays become severely limited 
by a lack of short baselines, and to cover astrophysically relevant dimensions between 0.1-3 milliarcseconds in 
visible wavelengths, baselines between pairs of telescopes should cover the whole interval 30-2000 m. 
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1. IMAGING STARS WITH INTENSITY INTERFEROMETRY 

1.1 Introduction 

A small number of the apparently largest stars have now been imaged with phase/amplitude interferometers, 
revealing their oblate shapes caused by rapid rotation or perhaps being surrounded by circumstellar winds, disks, 
or shells. Clearly, improving the angular resolution by just another order of magnitude would enable many more 
phenomena to be resolved. 1 However, atmospheric turbulence makes the required interferometry over hundreds 
of meters quite challenging for phase/amplitude interferometers which require optical path- length precisions 
within a fraction of an optical wavelength. 

Optical intensity interferometry currently seems the most realistic way to realize imaging of stars on scales 
well below one milliarcsecond (mas). 2 ' 3 Using a simple X/D criterion, 1 mas resolution at A=500 nm requires 
a baseline around 100 meters, while 1 km enables 100 /ias. For the forthcoming large arrays of Cherenkov 
telescopes, extensions over some 2 km are discussed, and if such could be utilized at A=350 nm, resolutions could 
approach 30 /ias. 

1.2 Theory of Intensity Interferometry 

Intensity interferometry measures correlations in light intensity fluctuations between pairs of telescopes. 4 ' 5 This 
correlation is related to the squared modulus of the complex degree of coherence of the light: 
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where 1{ and Ali are the intensity and intensity fluctuations respectively, measured by telescope i and () denotes 
the time average. 

The complex coherence is proportional to the Fourier transform of the stellar surface brightness, B, sampled 
at a point (u, v) in the Fourier plane, according to the van Cittert-Zernike theorem: 

7(u, v)<x J J 9(Z, m ) e - 2 ^ ul + v ^ dl dm (2) 

where /,m are (angular) positions in the sky. So, in short, each pair of telescopes in an intensity interferometer 
measures the square of the magnitude of the Fourier transform of the target at one point in the Fourier plane. 

Long telescope separations, or baselines, measure high-frequency components, corresponding to small struc- 
tures on the target, while short baselines sample the low frequencies. For a baseline B = (£?North, ^East) the 
associated Fourier-plane coordinates are (u,v) = ^(^North, £>East)- 

For an array of N telescopes there are N(N — l)/2 (not necessarily unique) baselines. As the target of 
observation moves across the sky, the projected baselines, B p , will change, with each telescope pair tracing out 
an ellipse in the Fourier plane according to the following expression: 6 

^ ^ / — sin / sin h cos h cos / sin h 
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where / is the latitude of the telescope array and 5 and h are the declination and hour angle of the star. The 
w component corresponds to the time delay between the two telescopes. This effect, known as Earth rotation 
synthesis, greatly increases the Fourier-plane coverage of a telescope array. 

From the Fourier magnitude, it is possible to both reconstruct the image of the target 7 (provided the sampling 
is good enough), or fit a model to the data. 
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Figure 1. Approximated angular diameters of the stars in the Bright Star Catalogue (BSC), 8 containing all stars in the 
sky with my < 6.5 (a few outliers have been omitted for clarity). Effective temperatures were estimated from their B — V 
color index, and then used to calculate the angular diameters. 



Optimizing the telescope placement in an interferometric array is a complex problem that has been widely 
studied for optical and radio interferometers (e.g., Refs. 9-11). The optimal configuration depends upon the 
task: if the objective is to study very small structures, long baselines are required. If, on the other hand, 
one wishes to measure large-scale structures, shorter baselines are necessary. For a general-purpose array, it is 
usually preferable to have a healthy mixture of short and long baselines, sampling the (u, v)-plane as uniformly 
as possible. In some interferometers this is satisfied by movable telescopes; in others one maintains similar 
diffraction-limited resolutions for observations at different wavelengths, placing telescopes in self-similar patterns 
(e.g., logarithmic spirals), while still others may be constrained by the local geography. In any case it is almost 
always best to avoid regularly spaced telescopes, since such placement results in duplicate baselines. 

1.3 Astrophysical Targets 

Ideal targets for intensity interferometry are hot and bright stars. 1 The limiting magnitude depends on the size 
of the flux collectors, the detectors with their electronics, and the integration time, but simulations and previous 
experience suggest a limiting magnitude around my ~ 9, 2 while initial observing programs will likely be limited 
to my < 6. 

The Bright Star Catalogue (BSC) 8 lists all stars brighter than my = 6.5 (a total of 9096), and can be used 
to illustrate the parameters of likely observational targets. The B — V colors of the BSC entries were converted 
to theoretical effective temperatures by fitting a polynomial to data relating T e ff and B — V color index from 
Bessel et al. 12 These effective temperatures were then used to calculate the angular diameters assuming the stars 
are blackbodies with uniform circular disks; see Fig. [I] 

Most bright stars have angular diameters on the order of one milliarcsecond (mas) or less. Such a star is 
resolved by telescopic baselines of ~100 m, so an array must provide baselines at least this long to be able to 
observe any significant number of objects. However, it is important to note that short baselines are also crucial. 
An array lacking baselines shorter than ~ 100 m will have too narrow a field of view to image the shapes of most 
stars. 

1.4 Technical Requirements 

Intensity interferometry places substantial technical requirements on the detection equipment used (for a dis- 
cussion of its technical aspects, see Ref. 13). Very high time-resolution detectors are required to measure the 
rapid intensity fluctuations, and fast hardware is needed for calculating correlation functions. However, by far 
the most expensive components are the flux collectors. The technique is sensitive to the square of the received 
light intensity, and large amounts of starlight need to be collected. On the other hand, there is no need for 
high-quality optics - a time resolution on the order of 1 ns, say, corresponds to a light travel distance of 30 cm, 
and optical imperfections comparable to such a magnitude can be tolerated. Despite the tolerance for optical 
imperfections, the cost associated with a large array of ~10 m telescopes is substantial, so a very attractive 
possibility is to operate intensity interferometry as an auxiliary program on arrays set up for other purposes. 

Seemingly ideal for this purpose are those air Cherenkov telescopes that are being erected for gamma-ray 
astronomy (e.g., Refs. 14,15). These measure the feeble Cherenkov light in air produced by cascades of secondary 
particles initiated by highly energetic gamma rays. They must have a time resolution of a few nanoseconds (set 
by the duration of the Cherenkov light flash); they must be sensitive to short optical wavelengths (Cherenkov 
light is bluish); they must be large (Cherenkov light is faint), and they need to be separated by distances of 
~50 m to ~200 m (extent of the Cherenkov light-pool onto the ground). 

2. EXISTING CHERENKOV TELESCOPE ARRAYS 

A number of air Cherenkov telescope facilities are already in operation, with the major ones summarized in 
Fig. [2] and Table (TJ 



Table 1. Properties of some existing Cherenkov arrays. N denotes the number of telescopes, A is the light collection area of 
each telescope, n is the number of unique baselines simultaneously available, B m in, B m ax indicates the range of baselines 
for observation at zenith. The corresponding range of angular diameters in milliarcseconds (1.22A/6) for observations at 
400 nm is indicated by m in,0max- 
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CANGAROO-III 16 
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57 
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100, 184 


0.5, 1.0 
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21 


50, 100 


1.0, 2.0 


H.E.S.S.-I 18 
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0.6, 0.8 


MAGIC 19 
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85, 85 


1.2, 1.2 


PACT 20 


25 


9.5 


300 


20, 128 


0.8, 5.0 


TACTIC 21 
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7 
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17, 20 


5.0, 6.0 


VERITAS 22 


4 


113 


6 


80, 140 


0.7, 1.2 



2.1 Interferometric capabilities of existing arrays 

Already the baselines offered by some of the existing arrays may be adequate for starting intensity interferometry 
experiments. Their adequately large flux collectors should provice excellent signal-to-noise ratios, but the longest 
baselines do not reach more than some 200 meters and their achievable angular resolution essentially overlaps with 
that feasible with existing phase interferometers (though observations could be made in the blue or violet, where 
the contrast of many stellar features is expected to be higher). Table [I] lists some of their relevant properties. 

Figure [3] shows a simulation of the Fourier plane coverage of the baselines possible with the VERITAS array 
for an 8-hour observation of a star moving through zenith, cf. Eq. While the (u, v) -cover age is increased 
greatly through Earth rotation, it is still rather sparse and does not extend very far. 
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Figure 2. Overview of existing air Cherenkov telescope arrays. The geometric layout is shown with north upwards (the 
scale varies between different arrays in the figure; see Table [l] for array sizes;. 

3. FUTURE CHERENKOV TELESCOPE ARRAYS 

The largest current projects for future observatories are CTA (Cherenkov Telescope Array), 23 and AGIS (Ad- 
vanced Gamma-ray Imaging System). 24 Either of them is envisioned to comprise around 50-100 flux collectors 
with diameters up to ~20 m spread over kilometer-scale distances. These would be excellent facilities for inten- 
sity interferometry, both in terms of resolution and Fourier-plane sampling, and their design studies now also 
incorporate their possible use for optical interferometrry. For such use, appropriate data analysis software would 
digitally synthesize baselines between all possible pairs of telescopes. 2,3 



100 




u-A (m) 

Figure 3. Simulated (u, v)-coverage of the VERITAS array (Fig.[2]and TableJT]) for an 8-hour observation of a star moving 
through zenith. Each pair of telescopes traces out two partial ellipses in the Fourier plane, since the Fourier magnitude 
in a point (u, v) will be the same as in the point (—u, —v). The filled circles indicate the snapshot Fourier plane coverage 
when the star is in the zenith. The small crosses mark 2-hour intervals before and after the meridian passage. 

3.1 Candidate configurations 

Any dedicated interferometer can be optimized for the best coverage of the interferometric (u, u) -plane. However, 
for Cherenkov telescope arrays, intensity interferometry is a secondary application and telescope placements 
and sizes will primarily be chosen for their gamma-ray observations. Nevertheless, also the interferometric 
performance may be optimized, without compromising the gamma-ray performance nor adding extra costs. 

A number of candidate array layouts for the CTA are being studied within its design study, 23,25,26 of which 
a number is shown in Fig. [4] Each of these candidate layouts is a subset of an all-encompassing larger array, 
shown in the bottom-right corner. Note that several among the configurations are not fundamentally different 
from each other. Already from a quick glance, one can roughly divide the layouts into three groups: A, B, F 
and G, whose telescopes are concentrated in a fairly small area and which lack longer baselines; C, D and J 
which contain telescopes more sparsely spread out over a large area; and E, H, I and K which contain a larger 
number of telescopes and provide both long and short baselines. For this paper, the configurations B, D and I 
were selected for more detailed study, as "representatives" of the three groups. 

The left column in Fig. [5] shows the telescope placement of the three examined configurations B, D and I. 

Table 2. Properties of the three examined candidate configurations. N is the number of telescopes, A is the light collection 
area of each type of telescope, n is the number of unique baselines available, B m in, Bmax indicates the range of baselines 
for observations in zenith. The corresponding range of angular diameters in milliarcseconds (1.22 X/b) for observations at 
400 nm is indicated by 6 m in, Omax- 
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Figure 4. Different telescope array layouts currently being discussed for CTA. Each of the configurations is a subset of the 
large array shown in the bottom-right corner. For this paper, configurations B, D and I were selected for further study. 

The middle column shows a "snapshot" (u, v)-plane coverage of the arrays for a star in the zenith. In the right 
column, the effect in the (u, v)-plane of the star moving 20 degrees from the zenith to the west is shown. 

Already from Fig.[5]and Table [2] something can be said about the designs if applied to intensity interferometry. 
Configuration B samples the central parts of the Fourier plane very densely and would provide a field of view 
that is fully capable of imaging the shapes of almost any stellar object. However, since it lacks longer baselines, 
it is unable to resolve details smaller than ~0.13 mas. 

Configuration D, on the other hand, provides baselines out to 2180 m, permitting studies of very detailed 
structures, down to around 50 /xas. However, its shortest baseline is 170 m, which means that any structures 
larger than ~0.6 mas will be lost. This makes it unsuitable for most stars (cf. Fig. [I]). Configuration I seems 
to provide the best of the two worlds: it has baselines short enough to measure the shapes of the disks of most 
stellar objects, while still providing very long baselines and a very good resolution. 



4. SIMULATED OBSERVATIONS WITH CHERENKOV TELESCOPE ARRAYS 

An intensity interferometer using two photon-counting detectors A and B and a digital correlator will measure 
the correlation of Eq. Q in a discrete form: 

Nab 



7(2) 



N A N B J 



-N, 



(4) 



where N A and Nb are the number of photons detected in A and B respectively, Nab is the number of joint 
detections and N is the number of sampled time intervals. In the simulations presented below, each measurement 
was simulated by generating random numbers N A , Nb and Nab an d inserting into Eq. Q . These will be Poisson 
distributed random variable^] with mean values = Pa • N, Xb — Pb * N and Xab = Pab * N. Here, Pa and 

*In reality, the measurement time is always long enough for the Poisson distributions to be adequately approximated 
as normal distributions. 
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Figure 5. Left column: Telescope placement for the three candidate configurations B, D and I (top to bottom). Middle 
column: (u, v) -plane coverage at an instant in time for a star in the zenith. Upper right-hand squares show the central 
400x400 m area. Right column: (u, v)-plane coverage for a star moving from the zenith through 20 degrees to the west. 



Pb are the probabilities of detecting a photon in A and B respectively within a small time interval At and Pab 
is the probability of a joint detection within At. 

These probabilites can be written out in terms of variables depending only on the instrumentation and the 
target of study (e.g. Ref. 27): 

P A = a A (lA)At (5) 
P B = a B (I B )At (6) 

Pab = PaPb + ^^{^{Ib^Iab? ^At 2 (7) 

Here a denotes the quantum emcienc}J^] of the detectors, (/) is the mean light intensity, r c is the coherence time 
of the light (determined by the wavelength and optical passband) and jab is the degree of optical coherence 
(proportional to the Fourier transform of the target). Such simulations were carried out for the various candidate 
configurations using the binary star in Fig. [6] as the source to be observed. 

4.1 Varying target magnitude 

Figure [7] shows the results of varying the brightness of the target. The simulations started with the star in zenith 
and then moving 50 degrees to the west in steps of 5 degrees, with an integration time of 2 hours at each time 
step (thus, such an observation would require several nights of integration). The observational wavelength was 
500 nm and the detectors were assumed to have a time resolution of 1 ns and a quantum efficiency of 70 %. 



The quantum efficiency of a detector is essentially the ratio of detected photons to incoming photons. 
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Figure 6. Test image of a binary star used for simulations and its (logarithmized) Fourier magnitude. 
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Figure 7. Simulated observations of the binary star in Fig. [6] with configuration B (top row), D (middle) and I (bottom). 
Vertical columns are for differently bright stars; left is for my = 3; center for my = 5, and my = 7 to the right. 



The data from the simulations are shown in the Fourier plane. To better illustrate the locations of low- 
amplitude structures and of noise, the data were logarithmized and normalized after first performing a linear 
interpolation between the "observed" points, in order to obtain the Fourier magnitude over a regular grid. A 3D 
shading was then applied to the data to further highlight the locations of structures (the absolute values of the 



measured correlations are not of any real significance in this context). 

These simulations are consistent with earlier results in that the practical limiting magnitude seems to lie 
somewhere around my = 6 — 7. Configuration B performs well in reproducing the three central lobes of Fig. 6(b) , 
but it lacks the longer baselines needed for higher- frequency components. 

Configuration D, on the other hand, measures the higher- frequency components very well but even for a 
target with such a small angular extent, the central parts of the Fourier plane - which determine the overall 
shape of the object - are distorted due to poor sampling at short baselines. Configuration I samples the whole 
Fourier plane well, but its smaller-size telescopes (see Fig. |5| makes it a bit more sensitive to target magnitude. 
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Figure 8. Subsets of the candidate configurations (B, D, I from top to bottom). In the leftmost column half the telescopes 
of the superset configuration were selected in a pseudo-random fashion. In the middle column, one in four telescopes was 
selected, and in the rightmost column one in eight. 



4.2 Observations with subsets of the configurations 

It is interesting to examine the performance of also various smaller subsets of the telescope arrays. It is not 
known whether some modifications or special equipment might be required at each telescope to be used for 
intensity interferometry, nor whether the full array will always be available. Observations using only a subset of 
the telescopes may thus represent a realistic mode of operation, at least initially. 

For each of the configurations B, D and I, three subset configurations were generated, as shown in Fig. [8j 
The configurations shown in the leftmost column, designated B2, D2 and 12, were obtained by selecting half of 
the full set of telescopes in a semi-random manner (attempting to preserve the overall "shape" of the array). 
In the middle column, one in four telescopes was retained (these are referred to as B4, D4, and 14) and in the 
rightmost column only one telescope in eight was kept (B8, D8, and 18). 



Figure [9] shows the output from simulations of the binary star in Fig. |6 (a)] using these subsets. The simulation 
parameters were the same as in Fig. [7] except that the magnitude of the star was now fixed to my =5 and the 
integration time was increased to 10 hours per timestep in order to depress measurement noise and thus highlight 
sampling effects for the various configurations. 
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Figure 9. Simulated observations of the binary star in Fig. [6] with subsets of configuration B (top row), D (middle row) 
and I (bottom row). The left column is for subsets containing half of the telescopes; the center for a quarter, and the 
rightmost for one eighth (see Fig.j8|. 

It is obvious that more telescopes are better in terms of Fourier plane sampling, and too few telescopes can 
ruin the results, as illustrated by the observation with D8. It should be noted, however, that no attempt was 
made to optimize the selection of telescopes in order to improve the sampling. Once the design for CTA is 
finalized, this could be one of the first problems to examine in order to determine which telescopes to equip first 
for intensity interferometry. 

4.3 Image reconstruction 

The discussion above has only dealt with Fourier magnitudes, which is what is directly measured by an intensity 
interferometer. Fourier magnitudes alone can be used to fit model parameters such as stellar diameters, binary 
separations, circumstellar disk thicknesses, etc. As is shown in Ref. 7, model-independent image reconstruction 
from the Fourier magnitude alone is indeed possible (though with some limitations). 



Figure 10 shows a proof-of-concept example of an image reconstructed from simulated data using one of 
the candidate configurations. However, image reconstruction is an active area of research and many practical 
issues remain to be improved. Thus, at this stage, it is not yet possible to tell which imperfections are due to 
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Figure 10. Image reconstructed from simulated data corresponding to 50 hours of observation of a magnitude 6 rotationally 
flattened star with an obscuring disk. The simulation was done using a telescope array qualitatively similar to those 
discussed in this paper (see Ref. 7 for details). The image in the bottom-right corner was used as input for the simulation. 
Also shown is the pristine image convolved with the point-spread-function of the array. The difference in sharpness 
between this and the pristine image is due to the inability of a finite-sized array to measure all the higher-frequency 
Fourier components. The convolved pristine image is thus the best image that is theoretically possible to reconstruct. 



fundamental limitations of the algorithms used or simply due to practical issues that can be solved. There is thus 
a risk that a comparison between reconstructed images from the different candidate arrays would be misleading, 
which is why the focus here is instead on the algorithm-independent Fourier magnitudes. 

5. CONCLUSIONS 

A facility such as the CTA will most likely be able to produce excellent scientific data if used as an intensity 
interferometer, providing measurements of spatial structures on scales down to a few tens of /ias - orders of 
magnitude better than what is possible today in visible light. 

The candidate configurations examined here - while not optimized for intensity interferometry - all provide 
very dense sampling of the (u, v)-plane due to the sheer number of telescopes. Apart from the simulations 
described above, a number of other tests were carried out to determine effects from varying the declination of the 
star and the geographic orientation of the telescope array. It was found that the very large number of telescopes 
make the effects of any such variations almost negligible. 

However, it must be noted that an array such as configuration D will be severely crippled by its lack of short 
baselines, essentially limiting such a facility to only studying stars smaller than ~0.5 mas. Also, to improve 
coverage of the (u, v) -plane and to ensure a wider variety of baselines (enabling better image reconstruction), 
an exact east-west orientation should be avoided for the grid of many smaller telescopes (since stars rise in the 
east, and move towards west, the projection of any east-west baseline will vary only in amount, not orientation). 
For the fewer large telescopes, repetitive geometric patterns should be avoided, in order to avoid redundancy in 
baselines. To cover astrophysically relevant dimensions between 0.1-3 milliarcseconds, baselines between pairs 
of telescopes should as far as possible cover the whole interval 30-2000 m. 

ACKNOWLEDGMENTS 



The work at Lund Observatory is supported by the Swedish Research Council and The Royal Physiographic 
Society in Lund. S. LeBohec acknowledges support from grants SGER #0808636 of the National Science Foun- 



dation. 



REFERENCES 

Dravins, D., Jensen, H., LeBohec, S., and Nunez, P. D., "Stellar intensity inter ferometry: Astrophysical 
targets for sub-milliarcsecond imaging," Proc. SPIE 7734, 7734-9 (2010). 

Le Bohec, S. and Holder, J., "Optical intensity interferometry with atmospheric Cerenkov telescope arrays," 
Astrophys. J. 649, 399-405 (2006). 

LeBohec, S., Barbieri, C, de Witt, W.-J., et al., "Toward a revival of stellar intensity interferometry," 
Proc. SPIE 7013, 70132E (2008). 

Hanbury Brown, R., [The Intensity Interferometer], Taylor & Francis, London (1974). 

Labeyrie, A., Lipson, S. G., and Nisenson, P., [An Introduction to Optical Stellar Interferometry], Cambridge 

Univ. Press, Cambridge (2006). 

Segransan, D., "Observability and UV coverage," New Astron. Rev. 51, 597-603 (2007). 

Nunez, P. D., LeBohec, S., Kieda, D., et al., "Stellar intensity interferometry: Imaging capabilities of air 

Cherenkov telescope arrays," Proc. SPIE 7734, 7734-47 (2010). 

Hoffleit, D. and Warren, Jr., W. H., "Bright Star Catalogue, 5th revised ed.," VizieR Online Data Cata- 
log 5050 (1995). 

Keto, E., "The shapes of cross-correlation interferometers," Astrophys. J. 475, 843-852 (1997). 

Millour, F., "All you ever wanted to know about optical long baseline stellar interferometry, but were too 

shy to ask your adviser," New Astron. Rev. 52, 177-185 (2008). 

Boone, F., "Interferometric array design: Optimizing the locations of the antenna pads," Astron. & Astro- 
phys. 377, 368-376 (2001). 

Bessell, M. S., Castelli, F., and Plez, B., "Model atmospheres broad-band colors, bolometric corrections and 
temperature calibrations for O - M stars," Astron. & Astrophys. 333, 231-250 (erratum 337, 321) (1998). 
LeBohec, S., Adams, B., Bond, I., et al., "Stellar intensity interferometry: Experimental steps toward 
long-baseline observations," Proc. SPIE 7734, 7734-48 (2010). 

Hinton, J., "Ground-based gamma-ray astronomy with Cherenkov telescopes," New J. of Phys. 11, 055005 
(2009). 

Volk, H. J. and Bernlohr, K., "Imaging very high energy gamma-ray telescopes," Exp. Astron. 25, 173-191 
(2009). 

Kubo, H., Asahara, A., Bicknell, G. V., et al., "Status of the CANGAROO-III project," New Astron. 
Rev. 48, 323-329 (2004). 

Chitnis, V. R., Acharya, B. S., Bhat, P. N., et al., "Status of HAGAR, the High Altitude Gamma Ray 

Observatory at Hanle," 29th International Cosmic Ray Conference Pune 5, 235-238 (2005). 

Hinton, J. A., "The status of the HESS project," New Astron. Rev. 48, 331-337 (2004). 

Cortina, J., "Status and first results of the Magic telescope," Astrophys. Space Sci. 297, 245-255 (2005). 

Bhat, P. N., Acharya, B. S., Chitnis, V. R., et al., "Pachmarhi array of Cherenkov telescopes," Bull. Astron. 

Soc. India 28, 455-457 (2000). 

Koul, R., Tickoo, A. K., Kaul, S. K., et al., "The TACTIC atmospheric Cherenkov imaging telescope," 
Nucl. Instr. Meth. Phys. Res. A 578, 548-564 (2007). 

Holder, J., Atkins, R. W., Badran, H. M., et al., "The first VERITAS telescope," Astropart. Phys. 25, 
391-401 (2006). 

"CTA, Cherenkov Telescope Array." http://www.cta-observatory.org (2010). 

"AGIS, Advanced Gamma-ray Imaging System." http://www.agis-observatory.org (2010). 

Bernlohr, K., Carmona, E., and Schweizer, T., "MC Simulation and layout studies for a future Cherenkov 

telescope array," Proc. 30th International Cosmic Ray Conference 3, 1469-1472 (2008). 

Hermann, G., "Towards the future Cherenkov telescope array CTA," Nucl. Instr. Meth. Phys. Res. A (2010). 

Mandel, L. and Wolf, E., [Optical Coherence and Quantum Optics], Cambridge University Press, Cambridge 

(1995). 



